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Phase diagrams of the one-dimensional extended Hubbard model with a nearest neighbor 
interaction V and site-off-diagonal (bond-charge X and bond-bond W) interactions at half- 
filling are obtained numerically, using the level-crossing ap proach for finite-siz e clusters. This 
analysis is straightforward extension of the preceding work |cond-mat /9909277 1 , where only the 
X term is considered as a site-off-diagonal term. Similarly to the X term, the W term enlarges 
bond-charge-density-wave (BCDW) and bond-spin-density-wave (BSDW) phases for W < 
and W > cases, respectively. In the strong-coupling region, however, the W term causes 
ferromagnetism for W > 0, and shifts the boundary of the phase separation in the U/t <C —I 
region, which are different from the roles of the X term. In addition, we discuss a parameter 
region where the BSDW state is the exact ground state, and show that the exactly solvable point 
appears just on the phase boundary between the BSDW and the ferromagnetic phases. 

KEYWORDS: site-off-diagonal interactions, g-ology, level-crossing approach, strong-coupling theory, bond-charge- 
density wave, bond-spin-density wave, exact solution 
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§1. Introduction 

The Hubbard model has been considered to be cm, 
the simplest models to describe electrons in metals.Erl 
However, since this model is derived by neglecting el- 
ements of the Coulomb integral except for the on-site 
term,El it may be unsufficient to describe phenomena 
where the neglected interactions play essential roles. The 
most dominant interaction next to the on-site term is 
considered to be a nearest-neighbor density-density in- 
teraction V. The Hubbard model including the V term 
is usually called "the extended Hubbard model" . In ad- 
dition, site-off-diagonal elements of the Coulomb integral 
give bond-charge X and bond-bond W interactions. Re- 
cently, effects of the site-off diagonal terms are pointed 
out in the real materials.o' In this paper, we consider 
the role of these site-off-diagonal interactions in one- 
dimensional (ID) cases, where non-perturbative treat- 
ment is possible. 

For this purpose, we consider the following generalized 
Hubbard model in the half-filled band, 

U = -t ^ + U ^2 n^Uil + V ^ n i n i+l 

i i i 

+W ^(5 M+1 ) 2 + X B iti+hs (m- s + n i+h - s ) 



X' Bj i j_|_i )S nj i _ s n,*-|-i ) _. s , 



where 



(1) 



(2) 



The three body term in eq. (1) (X' term) is not derived 
from the generalized Coulomb integral. This term is jus- 
tified, for example, as an effective interaction in the three 
band modellf In fact, the W term is interpreted as an 



effective interaction in electron-phonon systems.© 

Analysis of this model (X' = 0) in the half-filled 
band n was performed by Campbell et al. ten years 
agoE'tP They discussed the existence of charge-density- 
wave (CDW), "bond-order- wave (BOW)" and ferromag- 
netic (FM) phases, but they did not discuss all the phases 
in this model. Especially, the nature of the "BOW" 
phase was left ambiguous. 

The property of this model can be clarified by the level- 
crossing approach which is based on the conformal field 
theory and the renormalization group. This method en- 
ables us to determine the phase boundaries with high 
accuracy from numerical data of the finite-size clusters. 
Using this technique, one of the authors obtained the 
phase diagram of the model (1) for W = oBlid' In this 
paper, we apply this approach to the case W ^ 0, and 
clarify the roles of the W term. 

The most important purpose of this paper is to clarify 
the nature of the "BOW" phase. Japaridze first pointed 
out using the bosonization technique that two types of 
"BOW" phases may appear in half-filleji. ID pSystems by 
the effect of the site-off-diagonal terms .OL!3tij The two 
"BOW" phases are bond-charge-density-wave (BCDW) 
and bond-spin-density-wave (BSDW) phases. Their or- 
der parameters are given by 



O 



0bsdwOO = (-1) j X SjB 



(3) 
(4) 



Recently, Itoh et al. discovered that the BSDW state 
gives the exact ground state in-some parameter region of 
a generalized Hubbard model. x3 In this paper, by using 
the same argument, we discuss the exactly solvable cases 
of the present model (1), and clarify the relation to the 
numerical results. 
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§2. Site-ofF-diagonal terms 

To begin with, we consider the role of the site-off- 
diagonal terms by intuitive discussion. First, we con- 
sider roles of the bond-charge interaction (X, X') terms. 
Throughout this paper, we fix the parameters as X' — 
—2X for simplicity. In this case, the bond-charge inter- 
action terms are rewritten as 



U > -2V + <5.9, U > Sg. 



*I>lci+i, s +H.c.)K_ s 



H+l, 



5 ) 2 - (5) 



This term conserves SU(2)<g>SU(2)/Z(2) symmetry of the 
Hubbard model under the following canonical transfor- 
mation, 



(6) 



Hereafter, we call eq. (g) the X term. In a hopping 
process of an electron with spin s, the hopping term and 
the X term gives —t when the two sites have zero or 
two electrons with the opposite spin —s, whereas they 
give — t + X when the two sites have one electron with 
spin — s. Therefore, we can expect that electron pairs 
are stabilized when X < 0. 

Next, we consider roles of the bond-bond interaction 
(W) term. This term can be rewritten as 



Hw 



-4Wj2( S i-Si+i+Vi-Vi+i-t), (7) 



where Si is the usual spin operator, and r\ i denotes the 
77-pairing operator whose components are defined by 

vt = (- 1 T c iAl' T li={-l) i cncn, = -{ni-l). 

(8) 

Thus, the W term includes the spin exchange and the 
pair hopping. Since a pair hopping process gives an en- 
ergy gain 2W , one can expect that the dimer state is 
enhanced for W/t < 0. In addition, for W/t > 0, one 
can expect appearance of ferromagnetism for U/t ^> 1 
and a phase separation (ferromagnetism in the 77-pairing 
space) for U/t -C — 1. 

Apparently, the W term also conserve the SU(2) <g) 
SU(2)/Z(2) symmetry of the Hubbard model, so that the 
present model keeps this symmetry when V = 0. The 
SU(2) symmetry of the 77-pairing is broken when V =/= 0, 
and the entire symmetry is reduced to U(l)<g)SU(2). 

§3. Weak-coupling region 

The instability of the model in the weak-coupling limit 
can be discussed based on the g-ology.t3 The low-energy 
behavior of eq. (1) at half- filling is described by two 
sine-Gordon models for the charge and the spin sectors. 
Then, the phase boundaries are determined in terms of 
the bare coupling constants. According to the results, 
the spin sector with SU(2) symmetry has a gap for 



U < 2V - Sg, 



(9) 



where Sg = 4X/ir + 8W. On the other hand, the charge 
sector with U(l) symmetry has a gap in the following 



region 



V > 0, U < 



(10) 



(11) 



This region is given by two Berezinskii-Kosterlitz- 
Thouless (BKT) transition lines. Note that the BKT 
line given by V = reflects the SU(2) symmetry of the 
77-paring. In the charge-gap region, there is a valley of 
the gap, and the gap vanishes on the Gaussian transition 
line, 



U = 2V + Sg, V < 0. 



(12) 



In the g-ology analysis, both X and W appear only 
through Sg, so that they play similar roles in the weak- 
coupling cases. For Sg — 0, CDW and spin-density-wave 
(SDW) phases appear in the charge-gap region, while in 
the metallic region, two phases appear where singlet or 
triplet superconducting (SS, TS) correlation is dominant. 
The BCDW appears for Sg < while the BSDW appears 
for Sg > 0, between the CDW and the SDW phases. 

Although the result of the g-ology is valid only in the 
weak-coupling limit, the description by the sine-Gordon 
model does not break down in the finite-coupling region, 
due to the concept of the Tomonaga-Luttinger liquid. In 
this case, the instabilities described by the sine-Gordon 
model can be identified as level-crossings in the exci- 
tation spectra in finite-size systems. In this approach, 
logarithmic corrections are canceled on the transition 
points, and the dependence of the system size L ap- 
pears as 0(1/L 2 ). This correction originates from the 
deviation from the linearized dispersion relation of the 
Tomonaga-Luttinger model. Therefore, phase diagrams 
can be obtained with high accuracy from the numerical 
data of the finite-size clusters. The detail of the level- 
crossing approach is discussed in ref. [ITj. 

§4. Strong-coupling region 

The differences between the role of the X and the W 
terms are clarified in the strong-coupling cases. First, 
we consider the transition between the phase-separated 
(PS) state and the SS phase, and the one between the 
SS and the CDW states. In the U/t < -1 region, the 
model can be mapped onto the spin-1/2 XXZ spin chain 
using the 77-pairing operators defiped-in eq. @ and the 
second-order perturbation theory :E3Et 



Wxxz — J 



where 



J 



A(t-X) 2 



>h 



1 

1 
(13) 



4W, J A = J + 4V. 



(14) 



In the XXZ spin chain, A = 1 gives a BKT-type tran- 
sition between the XY (spin fluid) and the Neel states, 
and A — — 1 gives the first-order transition between the 
XY and the ferromagnetic phases. In the present case, 
the PS, the SS and the CDW states correspond to the 
ferromagnetic, the XY and the Neel states, respectively. 
Therefore, V = gives the SS-CDW transition, and the 
asymptotic boundary of the SS-PS transition is given by 

2(1 -£,ft 2 



V 



\u\ 



+ 2W, 



(15) 
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where £ = X/t. Note that this result is valid only for 
W/t < 0, because the SS-PS boundary do not intersect 
the SS-CDW line (V = 0) due to the property of the 
XX Z spin chain. 

Next, we calculate the energies of the CDW and the 
SDW states by extending van Dongen's perturbative ex- 
pansion^ Since the CDW state is described as a "Neel- 
ordered phase" (A > 1) in the ^-pairing space, the en- 
ergy is given by 

^™ = | + MW \ £XXZ (W/W - 1\) - 



(16) 



L 2 1 1 ' 17 (3v-l)U 

(1 - £) 2 [(36v 2 -5v- 1)(1 - £) 2 - 8(3v - l)v] t 4 
+ v(3v - l) 3 (4t> - l)U 3 ' 

where v = V/U. £xxz(A) is the energy density of the 



XXZ spin chain (|13|) with J = 1, in the Neel-ordered 
region (1 < A), 



£xxz(coshf 



OO 

n— — oo 



sinh 9 



-2\n\6 ' 



(17) 



where cosh^ = A. On the other hand, the perturbative 
expansion for the SDW state can be performed up to 
the fourth order, using the Bethe-ansatz result for the 
spin- 1/2 Heisenberg chain as 

"(l-O 2 * 2 



Esuw 



V - 



9C(3) 



- W 
(l-v)U 

(1-0 2 [2(1- 



4 In 2 



(18) 



if - 1 + v ] t 4 



(l-v) 3 U 3 

The expansion for the SDW state is valid only for W < Q. 

The energies for the PS and the FM states are not 
affected by the perturbation in the thermodynamic limit. 
In the PS state, the system separatea.into two domains: 
doubly occupied sites and a vacuum.H2P Then, the energy 
in the thermodynamic limit is given by 

E PS = ^~^L. (19) 
On the other hand, the energy of the FM state is given 

by 



Eym = VL. 



(20) 



Using the energies obtained above, we determine 
asymptotic forms of the four first-order transition lines 
by equating the energies as 

(21) 
(22) 
(23) 
(24) 

In the U, V — > oo limit, the CDW-SDW phase boundary 
is V-U/2 = 2|VF|(21n2-l) for W < 0, while the CDW- 
FM phase boundary is V - U/2 = 2W for W > 0. The 
FM-PS boundary is fixed at U = -2V. 

In the present strong-coupling theory, X appears in 
the every order of t 2n /U 2 "^ 1 in the perturbation series, 
while W appears only in the zeroth order. Therefore, 



Ecdw 


= EgD W 


Ecdw 


= EpM, 


Esdw 


= Ep S , 


Efm 


= E PS . 



the W term plays roles to shift the boundary lines of 
the phase separation, the CDW-SDW and the CDW-FM 
transitions. 

§5. Exactly solvable case 

We discuss the exact ground state of the BSDW state 
in the present model, using the argument in ref . [l4|, where 
the three body term (the X' term) is not included. The 
explicit wave functions for the BCDW and the BSDW 
states are given by products of "electron-hole dimers" 
as, 

L/2 

|BCDW) S = Y[[2i - 1, 2i],[2t - 1, 2i]_ s |0>, (25) 

i=l 
L/2 

|BSDW) S = \{[2i - 1, 2i],[2i, 2t + 1]_,|0), (26) 
t=i 



where [i,j] s = -^( C L + C J S )- These "electron-hole 
dimers" can be interpreted as electrons located on the 
bonds. 

Now, we consider a local operator hij S for spin-s which 
satisfies hij S [i, j] s \0) = 0, and hij s \ otherwise) > for the 
other configurations. For example, 



h a - 



n 



n ; 



1.0 _ r> 

n ijs ~~ D ijs 



-B 



hijs — Bij s + B l]s 



His ^js 2, 
1, 

2 



(27a) 
(27b) 
(27c) 
(27d) 



If the total Hamiltonian is written in the following form 
using the local operators as 



n 



/ , A J n i,i+l,-[ n i,i+l,l' 



Xj > 0, (28) 



then W|BSDW) S = 0. Therefore, the BSDW state is the 
exact ground state of this Hamiltonian.H The present 
model (1) can be written using the above local operators 
as follows 

7~L — E = [^2^(fri,i+i,f^i,i+i,4, + ^M+i.T^M+i.i^ 

i 

+m,i + i,t h i,i+i,l + Ki+r,l h ti+u)] . ( 29 ) 

where < A < 1, and the parameters and the ground- 
state energy of the model are 



U = 2V = 2(1 - X)t, 
X = -X'/2 = -Xt, 
W= (l + X)t/2, 
E = (1 - X)tL. 



(30a) 
(30b) 
(30c) 
(30d) 



The Hamiltonian (1) has SU(2) symmetry under global 
rotation in the spin space. In other words, the opera- 
tor for the total spin S" ot = J^i commutes with the 
Hamiltonian, so that (S' t + t )' l |BSDW} s with n < L/2 also 
gives the ground state. Therefore, the fully ferromag- 
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netic (FM) state (n = L/2), 

|FM) = f[ C t|0), (31) 

i=l 

also degenerates with the BSDW state. This large de- 
generacy is consistent with the result of the bosonization 
theory where a BSDW isjiescribed as bond-located spins 
with gapless excitations. I3't3 

In finite-size systems, there is no spontaneous breaking 
of the translational symmetry, so that the actual wave 
functions of the BSDW states are reconfigured by linear 
combinations of eq. (26) for s =| and j as 

|BSDW)± = |BSDW) T ± IBSDW)^ (32) 

These two states appear at two different wave numbers 
k = 0, 7r. When the parameters of the model shifts from 
the exactly solvable case toward the spin-gapless region, 
the k = mode remains the ground state, while the 
k = 7r mode gives the lowest excited state. 

We can also construct a model which has a BCDW 
state as an exact ground state, if the local operators in 
eq. ( p8[ ) is replaced as h^ i+1 ^ — > h"^ In this case, 
the obtained Hamiltonian includes interactions among 
three sites. 

§6. Numerical results 

Using the level-crossing approach as in the same way 
of ref. |l0|, we determine the phase diagrams of eq. (1) in 
the U-V plain for W ^ 0, X = X' = by analyzing the 
numerical data of the exact diagonalization for L = 12 
systems. We determine first-order transition lines for the 
FM and the PS states by comparing the ground-state 
energies with eqs. ( p0| ) and (|l9|). Then, we obtain the 
results as shown in Figs. l(d)-(f). For comparison, we 
also show the phase diagrams for W = 0, X = —X'/2 ^ 
in Figs. l(a)-(c) which are taken from ref. [l(j. 

For W < 0, the BCDW phase appears between CDW 
and SDW phases, while for W > 0, a BSDW phase ap- 
pears instead of the BCDW phase. This property is sim- 
ilar to the effects caused by the X term. The result 
is consistent with the estimation that a dimcr (BCDW) 
state is stabilized for X, W < as was discussed in Sec. 0. 
In ref. |8|, existence of the "BOW" state is discussed for 
X = X' = 0, W > 0, so that this "BOW" phase corre- 
sponds to the BSDW state in the present results. 

The differences between the roles of the X and the W 
terms are clarified in the strong-coupling region. One is 
a property of the PS boundary in the U/t <C — 1 region. 
For W < 0, the PS boundary shifts toward the negative- 

V side, and flows into (U,V) = (-oo,-2\W\). On the 
other hand, for W > 0, the boundary shifts toward the 
positive-V side, but it does not cross the V = line. 
Therefore, a direct CDW-PS transition takes place at 

V = 0. 

The other difference is appearance of a FM state for 
W > 0. For W/t = 1/5, a tetracritical point appears 
among PS, TS, SDW and FM states. As W is increased, 
the SDW region and the Gaussian line in the charge part 
are swallowed by the FM region, and a direct BSDW-FM 
transition takes place. 



For W/t = 1/2, an exactly solvable point (A = 0) 
appears at (U/t, V/t) — (2, 1). This point is just on the 
boundary between the BSDW and the FM phases as was 
discussed in Sec. ||. The exactly solvable point moves on 
the U = 2V line as A is increased. When W/t — —X/t — 
1 (A = 1), exactly solvable point appears at (U ft, V/t) = 
(0,0). In this case, CDW, PS, FM and BSDW phases 
degenerates on this point. The other regions of the phase 
diagram are occupied only by CDW, PS and FM states 
(see Fig. 2). 

In order to check the validity of the result in the strong- 
coupling region, we show in Fig. 1 (d), asymptotic forms 
of the PS boundaries given by eqs. (Jlf) and (23). H They 
well agree with the numerical result. In Fig. 3, we show 
the phase boundaries near the U = 2V line for L = 8-14 
systems, and asymptotic boundaries given by the per- 
turbation theory up to the second and the fourth or- 
ders. In Fig. 3 (a), the CDW-SDW transition line given 
by eq. (21) agrees with the Gaussian line in the charge 
part. Although the size dependence of the the spin-gap 
phase boundary is large, it also approaches to the same 
line in the strong-coupling regime as the system size is 
increased. In Figs. 3 (b) and (c), eq. (22) agrees with 
the CDW-FM transition lines. In Fig. 3 (b), the BCDW 
phase appears in the very narrow region, and a tetracrit- 
ical point appears among the BSDW, the BCDW, the 
SDW and the CDW phases. 

§7. Summary 

We have studied the phase diagram of the extended 
Hubbard chain including the bond-charge X and the 
bond-bond W interactions, using the level-crossing ap- 
proach. We have clarified the similarity of the roles of 
the X and the W terms in the weak-coupling case (ap- 
pearance of BCDW and BSDW phases), and difference 
in the strong-coupling case (appearance of a FM state, 
properties of the PS boundary). We have also shown the 
exact BSDW state appears just on the phase boundary 
between the BSDW and the ferromagnetic regions. 
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Fig. 1. Phase diagram of the ID extended Hubbard model with 
site-off-diagonal interactions determined by the data of the L = 
12 system at half-filling. [CDW (SDW): charge- (spin-) density 
wave, BCDW (BSDW): bond-CDW (bond-SDW), SS (TS): sin- 
glet (triplet) superconducting phase, PS: phase-separated state, 
FM: ferromagnetic phase] The asymptotic phase boundaries for 
the PS are given by eqs. ( fi5| ) and (23). (a)-(c) are taken from 
rcf. [uj for comparison. 
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Fig. 2. Phase diagram of the ID extended Hubbard model with 
site-off-diagonal interactions (W/t = —X/t = 1) determined by 
the data of the L = 12 system at half-filling. 
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Fig. 3. Phase boundaries [(G) Gaussian transition for the charge 
sector, (S) spin-gap transition] near the U = 2V line for (a) 
W/t = -1/5, (b) W/t = 1/5 and (c) W/t = 1/2, determined 
in L = 8-14 systems. The strong-coupling results [eqs. (21) and 
(22)] agrees with the Gaussian line of the charge sector in (a), 
and the CDW-FM boundaries in (b) and (c), respectively. 



